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1. Introduction 

The radiation dose received from high energy galactic cosmic rays (GCR) is a limiting 
factor in the design of long duration space flights and the building of lunar and martians 
habitats. It is of vital importance to have an accurate understanding of the interactions of 
GCR in order to assess the radiation environment that astronauts will be exposed to. 

Most previous studies have concentrated on strong interaction process in GCR. However 
there are also very large effects due to electromagnetic (EM) interactions. EM studies have 
previously concentrated on single photon exchange leading to nucleon removal. However two- 
photon processes also occur which lead to the production of lepton pairs with cross sections 
of the order of kilobams. Also at high energy the stopping powers from these processes can 
exceed that due to atomic collisions. Thus even though very high energy GCR are not as 
abundant as lower energy GCR they still must be considered due to the fact that the cross 
sections and stopping powers are so much larger than normal. 

In this report we describe our first efforts at understanding these EM production processes 
due to two-photo collisions. More specifically, we shall consider particle production processes 
in relativistic heavy ion collisions (RHICs) through two-photon exchange. Examples of this 
broad category of processes include: 


z x z 2 -» ZiZ 2 /+r 

(1.1a) 

Z\Z 2 —* Z\Z 2 s + s~ 

(1.1b) 

Z1Z2 — ZiZ 2 V + V~ 

(1.1c) 

Z\Z 2 - ZiZ 2 tf° 

(lid) 


in which l + l denote charged leptons, s + s denote charged scalars, V + V~ denote charged 
vector particles, and is a neutral Hi gg s scalar. 

We shall limit our consideration to cases in which the colliding nuclei are identical, so 


that Z\ — Z 2 = Z. An important Feynman diagram that contributes to (1.1a), (1.1b), and 
(1.1c) is shown in the following figure (fig. 1). 
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For process (l.ld), an important diagram is shown in figure 2, in which the t riangular 
loop receives contributions from quarks, leptons and W gauge bosons. These processes are 
important for the following reasons (ref. 1). 

(1) These kind of processes become increasingly important as energy of the colliding nuclei 
increases, since their cross sections increase with energy. Thus their contributions to the 
stopping power of high energy ions also become more important at high energies. 

(2) These processes can be channels for production of charged particles, e.g., l + l~, W + W~, 
and neutral particles such as Higgs bosons, and various mesons. 

(3) For high Z nuclei these processes can be used for studying non-perturbative effects in the 
electromagnetic interaction. 

(4) They must be taken into account in the study of strong interaction effects in heavy ion 
collisions since they can lead to important background events, and must be taken into 
account also in the design of experimental set up, since they can lead to significant beam 
loss. 

Section 2 of this report gives a brief survey of a few major approaches used in the 
calculations for these processes. Section 3 examines some results of our calculations. We 
then point out briefly some open questions and make a few concluding remarks in Section 4. 

The purpose of this report is threefold. (1) It gives a simple, elementary introduction to 
this field. (2) It provides sample calculations for illustrating the approach we use. (3) The 
background and tec hni ques developed here can be used as a general base for la unchi ng further 
and more specialized studies into this field. 

While it is not our main goal here to obtain new and original results, some of our results 
are possibly new, and are as yet not available in the literature. 


2. A Brief Survey of Different Approaches 

In this section, we briefly list a few major approaches used in calculating cross sections 
for the kind of processes we are interested in. The first approach has been discussed in 
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references 2 and 3. In this approach, each colliding nucleus is replaced by an equivalent 
spectrum of photons. Each nucleus is considered to move in a straight line, unperturbed by 
the interaction. At a distance 6 from the line of motion of a nucleus, a spectrum of photons 
is generated, whose frequency distribution has the form: 




I (*) + ~2 K 0( x ) 


(2.1a) 


where 


uib 

jv 


(2.1b) 


Kq , K\ are modified Bessel functions, see reference 4, Sections 3.7 and 15.4. 

The cross section for this process can be written as an integral of a photon distribution 
function multiplied by a photon-photon cross section. 


= f— f —F(ui,u}2)a ry {ui,u 2 ), 

J w\ J 


(2.2a) 


where 


r oo roo y*2ir 

F(ui,u>2) = 27t / b\dbi I &2^2 / d0JV(u;i,6i) 

— ^1 mn 


'bl mn * £>2mn 

x -R1-R2) 


(2.2b) 


and 

b' = + b\ - 26162 cos <t >) 1/2 (2.2c) 

where uj\ and u 2 are the frequencies of the photons emitted by the nuclei, 61 and 62 are the 
distances of the nuclei from the point where the photons collide. Details can be found in 
Appendix A. Various differential cross sections can be derived from these equations. First 
we consider where W is the mass of the produced charged particle pair. We note that 
W 2 = 4m\<uj 2 . Hence we can equate in (2.2a) 


(2.2a) 


f dun fdw 2 / w 2 \ ( w 2 \ 

° = J -£■ J wt f r- * m 
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and 


da 
dW 2 


r) ffTT ( W1 ’^)' (2 ' 3b) 

Next we define the probability for producing a particle pair P(6) at impact parameter b by 

1 da 


p(b) = 


27 xbdb 


(2.4a) 


where 


^ = /— j (2.4b) 

ab J u\ J u>2 

in which it is understood that the 6 function is to be taken inside the triple integral which 
defines F(u/i,u; 2 ). The correctness of (2.4) can be checked by integrating both sides of (2.4b) 
over all values of the impact parameter b, which then yields (2.3a) for the total cross section. 
P(b) is the probability for the events in which two nuclei collide with each other at impact 
parameter 6, producing a charged particle pair in the process. A quantity L, known as the 
two-photon luminosity function is defined by (see ref. 2, eqs. (1), (9), and (10)) 

r2 / m/2 ' 


_ [dm [dw 2 ( 

J m J w 2F \ 1 ’ 4u>i ) 


So 


and 


dL 

dW 2 


1 f ^ i? ( ^ 2 \ 


da dL /rir2\ 

dw 2 ~ dw 2<Jrr(w 


(2.5a) 


(2.5b) 


(2.5c) 


where we have used the fact that ayy(u>it^2) actually depends only on W 2 so that we can 
write 

a ry (u 1 ,u 2 ) = a^iW 2 ). (2.6) 

It is our view that equation (10) of reference 2 is in error, and have duly corrected the error 
in the above definition of the luminosity function L. For stopping power calculation, we use 
the formula 

~~T~ =P f ~ /— (^1 +U2)F(u}i,U2)a ry {ui,U2), (2.7) 

dx J u>\ J ll >2 
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where p is the number of nuclei per unit volume. 

The second type of approach has been applied to a related set of purely quantum 
electrodynamic (QED) processes: e+e~ — * e+e - /"*"/ . This process can be calculated 
within the framework of QED. Cross sections can be obtained numerically by Monte-Carlo 
integration. Approximate formulas for total cross sections have also been obtained. See 
references 5 and 6. This kind of approach can be modified to apply to RHIC processes, 
provided one t ak es into account properly the effects of nuclear currents. See reference 7, 
Section II. 

In an approach closely related to this second type of approaches, Bottcher treated the 
colliding nuclei class i cally , by regarding them as classical charge distributions. The remaining 
amplitude for the production of charged particle pair is then obtained in the framework of 
QED. Thus for the case of the reaction Z\Z 2 -*• ZiZ 2 l+l~ the total cross section can be 
written in the form (ref. 8, eq. (10), p. 38): 


zjzj 

4t/ 2 


(A [ d p ~ 

(4 '°) J IS 


d3p - d3p+d%1± /{(kfjf 2 2 {4) 

h 


) 8 2p_o2p+o 

1 


+ A 


i>-~ fa - 


i>~- fa - m t 

2 

I 

«(p+, »+) 


( 2 . 8 ) 


where v denotes the velocity of one of the nuclei in the center of momentum frame, p— and 
p+ are the momenta of the produced leptons, s- and s+ are their polarizations, u(p-, S-) 
and v(p+, s+) are the lepton spinors, fci and fa are the momenta of the exchanged photons, 
and /i and f 2 are the nuclear form factors. For any 4-vector A, the slash notation 4 is 
defined by 

4=E Vr". < 2 - 9 ) 

/i=0 

where 7^ are the Dirac 7-matrices, see for instance reference 9, Appendix 2, pages 355- 361. 
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3. Results 


In this report, we adopt the approach discussed in references 2 and 3. As samples of our 
calculations, we present a number of results for the process 208 P6 208 P6 — ► 208 P6 208 P6 l + l~, 
and some others. Most of our calculations are done for colliding beam energies of 3400 Gev 
and 8000 Gev per nucleon. The impact parameter b varies over the range from 10 fm to 
1000 fm. The mass of I + l - varies from a threshold equal to 2m/ up to about 1000 Gev. In 
Appendix B we list the photon-photon cross sections for the following processes: 


77 - i+r 

(3.1a) 

1 

T 

(3.1b) 

77 -4 V + V~ 

(3.1c) 

ry~H° 

(3. Id) 

ss sections are also given there. 

By using (2.2)- 


(2.6), we can then obtain various luminosity functions, differential and total cross sections, 
probabilities, and stopping powers. 

Table 1 shows the total cross section for 2 ® 8 P6 2 ® 8 P6 — > 208 P6 208 P6 e + e - . We compare 
our numerical results based on (2.1) and (2.2), with the results based on the approximate 
formula (ref. 5, eq. (F.l), p. 276) 


a = 


- (l 3 - Al 2 + BI + C) 

27nm£ 


(3.2) 


where 


A « 6.36, B w 15.7, C ~ —13.8, 

Z = i n ('H£L_E?), 

V m 1 m 2 ) 

Zi, Pi, and mi, i = 1,2 are the changes, momenta, and masses of the colliding nuclei. 
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Table 1 


Incident energy /nucleon 
E(Gev) 

Total cross section (£m 2 ) 

Our results 

Calculated from formula 

3400.0 

8000.0 

0.8382 x 10 8 



Table 2 shows the corresponding stopping power calculations. The energies of the incident 
particles are given for both the case of colliding beams and also the case of an incident beam 
colliding with a fixed target. 


Table 2 


Incident energy /nucleon 
E (Gev, colliding beams) 

Incident energy /nucleon 
E (Gev, fixed target) 

KBKpH 

0.9636 

1.039 

0.2496 x 10 1 

0.1367 x 10 1 

3.039 

0.1918 x 10 _1 

0.2704 x 10 1 

14.64 

0.1725 x 10 -1 

3400.0 

0.2462 x 10 8 

0.8855 x 10 1 

8000.0 

0.1363 x 10 9 

0.3131 x 10 2 

For Pb-208, rho = 


In figure 3, we give plots of as a function of W in different ranges of W . The 

differential cross section can be obtained from by multiplying by a 77 cross 
section as in (2.5c). 

Figures 4a-d show plots of oy^W 2 ) for the reactions 77 — ► l + l~, 77 — > s + s“, 
77 — V + V~ , and 77 — H°. 

Figure 5 shows plots of P(b) for the reaction 208 Pb^Pb — *■ ^Pb^Pb e + e~ at different 
energies. 

Figure 6 presents plots of the total cross section for the process 208 P 6 208 P 6 — * 

208 Pb 208 Pb H°. 
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We have compared some of our results with the published results of Papageorgiu and Baur, 
and found some good agreement. In the following, we give a sample of such comparisons. 

Table 3 



W (Gev) 

W/J i 

W 2 .-#- 

dW? 

Our result 

Papageorgiu ’s 
result 

Incident 

energy/ 

nucleon 

(colliding 

beams) 

E = 3400.0 
Gev 

100.0 

0.7070 x 10 -4 

0.3152 x 10 3 

0.33 x 10 3 

141.4 

0.1000 x 10 -4 

0.8630 x 10 2 

0.90 x 10 2 

212.2 


0.1206 x 10 2 

0.13 x 10 2 

282.9 

0.2000 x 10 -3 

0.1990 x 10 1 

0.21 x 10 1 

E = 8000.0 
Gev 

19.0 

l 

0.5709 x 10" 4 

0.6129 x 10 3 

0.70 x 10 3 

280.0 

0.8413 x 10” 4 

[ 0.1708 x 10 3 

0.20 x 10 3 

370.0 

0.1112 x 10“ 3 

0.5444 x 10 2 

0.60 x 10 2 

460.0 

0.1382 x 10~ 3 

0.1881 x 10 2 

0.20 x 10 2 

550.0 

0.1653 x 10 -3 

0.6866 x 10 1 

0.70 x 10 1 

640.0 

0.1923 x 10 -3 

0.2606 x 10 1 

0.28 x 10 1 









































Table 4 


Incident energy /nucleon (colliding beams) E = 3755.6 Gev 

W(Gev) 

^7 (frn 2 Gev 2 ) 

Our result 

Baur’s result 

0.1200 xlO -2 

0.6188 xlO 13 

0.62 xlO 13 

0.1414 xlO -2 

0.4923 xlO 13 

0.48 xlO 13 

0.1732 xlO -2 

0.2777 xlO 13 

0.28 xlO 13 

0.2000 xlO -2 

0.1714 xlO 13 

0.17 xlO 13 

0.2200 x 10" 2 

0.1222 xlO 13 

0.12 xlO 13 

0.2400 x 10 -2 

0.8894 xlO 12 

0.92 xlO 12 

0.2600 x 10 -2 

0.6607 xlO 12 

0.69 xlO 12 

0.3000 x 10 -2 

0.3846 xlO 12 

0.40 xlO 12 


Papageorgiu and Baur’s results were taken from appropriate graphs in their papers (ref. 
2, fig. 3; and ref. 3, fig. 9). 

Cross sections are expected to scale roughly as Z 2 Z 2 . For our case Z\ = Z% = Z . So in 

order to obtain the corresponding cross sections, luminosity function, or stopping power for 

Z\ z 2 

different nuclei, one can simply multiply the results we have here by a factor (SV- Ttos 
if one wants the results for Al Fe collision, one can multiply the results presented in this 
section by the conversion factor The different nuclear sizes are expected to affect 

the results also. However for a rough order of magnitude estimate, such a simple scaling is 
expected to be reasonably accurate. 

4. Open Questions and Conclusions 

For small values of 6, and m/, such as mj = m e , P(b) exceeds 1. This signifies the 
breakdown of perturbation theory. The question as to how to extract meaningful results 
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from theory is under active investigation. See reference 10. In our simple approach, 
we have regarded the nuclei as point charges. By using form factors for the nuclei, the 
problem of violation of unitarity is expected to be somewhat ameliorated. However this 
problem still needs to be addressed, because for high Z nuclei, the coupling constant for 
the electromagnetic interaction is of the order Ze, even with nuclear form factor taken into 
account, which may therefore still lead to a breakdown of the perturbative approach to cross- 
section calculation. In a collaboration with Mirek Fatyka of Brookhaven National Laboratory 
(BNL), we shall investigate production and neutral meson production (such as c°, rp) in 
high energy heavy ion collisions. In these processes, we shall look for possible deviation in 
the measured rates or cross sections for values calculated by perturbation theory. 

In many studies of the type of processes considered here, various approximations are 
used. We have mentioned the equivalent photon approximation, and the semi-classical 
approximation. Also, in the approach of references 2 and 3, which we have adopted in 
this report, the effect due to phase coherence of the electromagnetic field generated by the 
nuclei has not been properly taken into account. One needs to investigate how valid these 
approximations are and what the regions of validity are for them. 

When one is primarily interested in the kind of electromagnetic processes disc uss d 
here, one needs to be able to estimate reliably the background due to strong interaction. 
Furthermore, there are other electromagnetic processes that also need to be studied, in 
addition to the ones we have looked at, even though the ones we have considered are among 
the most important. 

In su mm ary, we have given a brief introduction to two-photon exchange processes in 
high energy heavy ion collisions. Our calculations are based on an approach discussed in 
references 2 and 3. In view of the significance of this class of processes, and the many open 
questions that remain to be answered, we believe that further study in these areas will be 
valuable, not only for gaining a better understanding into these processes themselves, but 
also for studies and experiments in strong interaction physics. 
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In the following Appendices, we discuss the derivation of some of the formulas we have 
used. We look at the equivalent photon approximation in Appendix A and show how this is 
applied to the two-photon exchange processes in RHICs. Then in Appendix B, derivatives are 
given for some 77 cross sections. Appendix C provides a derivation of fermion contribution 
to the process H° 77. In Appendix D, we look at the details of how certain integrals 
encountered in our calculations are evaluated. Finally Appendix E gives a simple derivation 
of the formula (2.7) used for calculating stopping power. 


Appendix A. Equivalent Photon Approximation 

Consider a charge of moving along the x-axis. The effect of this charge on another charge 
located a distance b from the x-axis can be approximately calculated as follows. 

By first considering the electromagnetic (EM) field due to q in its own rest frame, and 
then making a Lorentz transformation to the laboratory frame, it is straight forward to show 


that the electromagnetic field due to q is given by 

Ei = -qv'/t(b 2 + 7 2 vV )~ 3 / 2 (A. la) 

E 2 = qlry(b 2 + 7 Vt 2 )~ 3 / 2 (A. lb) 

£3 = ^E 2 = q- 67 (ft 2 + 7 2 v 2 t 2 ) -3 / 2 (A.lc) 

c c 

E 3 = Bi=B 2 = Q (A. Id) 


t = 0 corresponds to the instant when q passes through the origin. When v « c, the 
components E 2 and £3 can be thought of as the components of a pulse of plane-polarized 
EM wave travelling along x. The energy flux of this EM field is given by the Boynting vector 

s= £ s * s - < A - 2 > 

So ignoring E\ for the moment, S points along x, and its magnitude is 

1^1 = 4^ E h (A.3) 

in which we have made the approximation ^ ^ 1. Over a unit area, the flow of energy is 

/I 151 * = T* C ^ dt (A ' 4a) 

Using Parseval’s theorem, we therefore have 

J |.S| dt = — J \E2{u)^ du /, (A. 4b) 

where £2 is the Fourier transform (FT) of E 2l defined by 

3 -k II m c ~" * ' A - 5 > 
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Hence the quantity S^w), defined by 


5 2 («) = £ IB 2 MI 2 , (A.6) 

can be thought of as the energy per unit frequency per unit area of the EM field at frequency 
<jj generated by the moving charge q. To obtain the photon number per unit frequency per 
unit area at frequency w, we set n 2 (u/) = ^jS 2 (ui), since each photon has energy hu. For 
the function N 2 (u;), the dependence on the distance b is implicit. To make the dependence 
on 6 explicit, we can write instead 

n 2 (u;,&) = (A. 7) 


From (A. lb), we obtain 


Hence 



(A.8) 


(A-9) 


The remaining component E\, of the EM field can be complemented by a magnetic field so 
that they can be considered to form a pulse of plane polarized EM wave. The same treatment 
can be applied to these components, so that the energy spectrum can be similarly obtained 
as before. The result is 


SiM = 


C . - . . .9 C 1 1 


/ 9\ 2 2 / w 6\ 2 

\bJ ir J 


x Kq 



(A.10) 


The effect of this pulse is roughly ^ that of the first pulse. So at high velocity, the second 
pulse can be neglected when compared with the first pulse. 

In conventional treatment, the two pulses are then simply added together, so that the 
effect due to the original moving charge q is replaced by a spectrum of photons whose number 
density is simply the sum of the number densities from the two pulses discussed above. Thus 
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one set 


■ is (D'l (^) 


_2 

After identifying q = Ze, ^ = a, and noting 


^(-u/) = 5j(u;) for i = 1,2, 


(AH) 


(A. 12) 


the photon energy spectrum 


N{uj , 6) = S\ (w) 4 - ^(uj) 4- Si(-uj) 4- <?2 ( — w) 


Z 2 e 2 1 


= n 



0’fe)’Kfe‘) 

©•&)>(£) 



(A.13) 


Application of Equivalent Photon Approximation to Two-Photon 
Exchange Processes 

When two nuclei Z\ and Z2 collide with each other, their EM interactions can be studied 
in terms of the EM interaction of the spectra of photons emitted by the nuclei. The situation 
can be pictured as in figure A.2. 

The two photons 71 and 72 are considered as colliding head-on with each other. Taking 
a cross-sectional view perpendicular to the direction of motion of the nuclei, the situation 
can be pictured as shown in figure A. 3. 

From our previous discussion, the number of photons emitted by Z\ at P, whose 
frequencies are between u/j and u\ + dui 1, is n(ui,6i)duii b\ db\ d<f> 1, where n(u/i • hj) 
is defined by (A. 11). Similarly, the number of photons incident at P emitted by Z2 is 
71(0)2,62)^2 &2 ^2* Therefore the EM cross section for the collision of Z\ and Z2 

through two-photon exchange can be written as 


a 


J n(u/i,&i) n{u2,b2) (011,0)2)61 d<f>\ 62 db^ dfa 

x 6 {b — R\ — R.2)du)\ du>2 


(A. 14) 
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in which R\ and R 2 stand for the nuclear radii of Z\ and Z 2 , and the ^-function takes into 
account that when b < R 1 + R 2 , the two nuclei overlap, and the EM interaction is swamped 
by the strong interaction of the nuclei, and so one needs to restrict 6 to values > R\ + R 2 if 
one wants to look only at EM interaction. 

Since 6 = (6j 4- 6j - 26162 006 the integration / d<f>\ dfa in (A. 14) can be simplified 
if one integrates over <p\ and converts the integration over into an integration over <p: 


J dtp 1 d4>2 


—* 2 7T 



(A. 15) 


So (A. 14) can be rewritten as 


= 27t Jn M ) ”(^2,62) a r y{ui,U2)d(b - R\ - R 2 ) 


(A- 16) 


x 61 6di 62 db 2 d<p cLji du 2 - 

If one now substitutes for n(ui, b{), i = 1, 2, using (A.ll), one obtains (2.2). 

Concerning the cutoff for 61 and 62, we observe the following. (A. 14) involves an 
approximation, which consists of replacing the virtual photons emitted by Z\ and Z 2 with 
real photons 71 and 72. This approximation is valid only if the masses of the virtual photons 
Ai and A2 are small compared to the mass of the produced system W. (See ref. 4, Sections 6.1 
and 6.7). By the uncertainty relations, A* = t = 1, 2. Hence in order for the approximation 
in (A. 13) to be valid, we must have A^ < W, or j- <H^. Therefore, 



(A. 17) 


If 6j does not satisfy (A. 17), contribution to the cross section is small, and is generally 
considered negligible. See reference 5, Sections 6.1 and 6.2, and reference 11, Sections 7, 
7. 1-7.3. Another consideration for the values of 6j is that since we are interested in the 
effects of each nucleus acting as a single entity rather than as a collection of nucleons acting 
independently of each other, i.e., we are interested in the coherent effects of the collection of 
nucleons, we need to restrict 

bi > Ri. (A. 18) 


15 


So for reactions in which the Compton wavelength of the produced system is smaller than 
the nuclear radii, i.e., ^ < R{, we can set the minimum of 6* by 

bi mn = Ri- (A.19) 

This is the case for n~ and r + r“ pair function. But for e + e - pair production, the 
Compton wavelength of an electron ^ is >-Rj. So we set the minimum of by 

(A. 20) 



16 



Appendix B 

First we list the cross sections for the processes in (3.1): 77 — * l + l ~, 77 5+3 , 77 — ► 

V + V~, and 77 — H°. (See ref. 2 , eqs. (14)-(17), pp. 159, 160; and ref. 12 , eqs. (10), ( 11 ), 


p. 95.) 


< 7(77 




w 2 

- (1 + vi)>fi -y/] . 


(B.la) 


< 7(77 — » s + s^) = 


27 rcr 
~W 2 


(1 + V,) - (1 - 5 v,) In + V^pi)] . 

(Bib) 


< 7(77 - v+V 1 ) = 


8na 2 


^1 + —ty + 3t^ A — 3tv(l — 2 ty) In > (B.lc) 


where 


__ 4 m 2 
W 2’ 


4m‘ 


M r i/2 • * 


(B.2a) 


tv — ^,2 > ^ — a/ 1 4 t w . 

W is the total energy of the two photons in the center of momentum frame. 

ff (rr - n°) = ((W 2 - mff) 

where T can be written as (ref. 12, eq. (10), p. 95) 

Cp’GpTTVjj 2 


r = 




87t 3 \/2 

and I in turn has the form (ref. 12 , eqn. (11), p. 95) 

I = ^2 0 2 h + ^ 

9 l 


(B.2b) 


(B.3) 


(B.4) 


(B.5a) 
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Iq — Z [2Ag + Ag(4Ag — l)/(Ag)j , 


(B. 5 b) 


where for 


and for 


h = 2A/ + Aj(4A/ - l)/(Aj), 

(B. 5 c) 

I w = 3A W r(l - 2 X w )f(X w ) - 3A^ - 

(B. 5 d) 

A>i, /(A) = -2 (arc^n^) , 

(B.6a) 

<>> = K ta £) -T + iTln $' 

(B.6b) 

” ± = \ ± 'J\- X 

(B.6c) 


The subscripts q, l, and W stand for quark, lepton and W-boson, respectively. 


m 2 

Xi = — for i = q,l,W, 
mi 


(B.7) 


and 77ij are the rest masses of the coresponding particles, is the rest mass of H G . 
In the following, we give the derivatives of the cross sections for the processes 


77 — s + s , 


(B.8) 


77 — I+r. (B. 9 ) 

We also give a derivation of the relationship between a and T for the process 

77 — H° (B. 10 ) 


77 — s + s . 
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The lagrangian for the system, including EM interaction, can be written as 


£ cm = - + ieA^j <f> + - ieA ^ <f> - m 2 0 + 0 (B.ll) 

in which <f> denotes a scalar field operator, Ay, denotes the photon field, p = 0, 1,2, 3. We 

use the convention that repeated indices are summed over, so that for example, 

Ay A* = AoA° + AiA 1 + A 2 A 2 + A$A Z (B.12a) 

= AoAq — A\Ai — A 2 A 2 — A$A$. (B.12b) 


This lagrangian can be separated into a free part, and an interaction part, so that 

0 


Cml - ie A<1 ^ ~ 


(B.13) 


+ e 2 A»Ay4>+4>. 

The 5-matrix element that contributes to (B.8) can be written in the form 


< p—;p+ |5^ + 5^| k\,ei;k 2 ,e 2 > (B.14) 

in which p± denotes the momenta of 5^, fci,€j are the momenta and polarization vectors of 
the photons, i = 1,2, and 5 ( 2 ) is defined by 

# 2 

= ~2 T J M <fix 2 , (B.15) 

where T denotes the time-ordering operator. Contribution from can be represented by 
the diagrams. 

Using (B.13) and (B.15), and standard techniques of field theory, one obtains 


<P-\ P+|£^|fcliCi;fc2>«2 > 

= te 2 (27r) 4 (2p_ 0 • 2 p+ Q • 2fci 0 • 2*2o) _2 ft -2 

+ P-/i) c 2 (p+^ ~ fey) , ^2 (M* +P~m) € i(P+^ ~ Kf) 

k 2 — m 2 + ie Jb' 2 — m 2 + ie j 

x 6*{p- + p+ - k\ - fc 2 )> (B.16) 
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where fc = p- — fci = fc 2 — p+, k' = p- — k 2 = k\ = p+, fi is the normalization volume, and 
e denotes an infinitessimal quantity. 

Likewise is defined by 


S (1) =*T J £ ml (i) J'k, 


(B.17) 


and 


<p_;p + |5 (1) |fci,ei;fc 2 ,C2 > = »e 2 (27r) 4 (2p_ 0 • 2p +a ■ 2k 10 • 2fc 20 ) 1/2 
x n -2 2ei e 2 6*(p- +p + -k\ -k 2 ). (B.18) 

The diagram representing this matrix element is shown in figure B.2. 

The total cross section is obtained by squaring (B.14), averaging over photon polarizations 
«1 and e 2 , integrating over phase space, and finally dividing by the photon flux. Hence we 
have 

|2 dPp~d?p+ 




<p-;p+ 


5(1) + 5(2) 


ki.ei; ^2>e 2 > 


(2tt)6 


„ 2 1 
xii z x — x 


(B.19) 


«lt«2 

r 

2c ~ T 0 ' 

in which T 0 is the normalization time. Substituting (B.16) and (B.18) into (B.19), we obtain 

' “ eV / /w 44 <"- + * s? 


= (P -” ~ mJ) ‘ /2 ^ ^ 4 M * 2lr ' < B ' 20 > 


where 


m = 




(A + p-)c 2 • (p +_ A) 


+ 


<i.«j ‘ * 2 -m 2 + « 

c 2 -(k / +p-)ei(p+ ~kQ 
k^-rr^ + ie 


i2 


+ 2ej • e 2 


a = (2 tt) 4 (2 p-o • 2p +0 ■ 2k w • 2*^) - 1 / 2 Q" 2 , 


(B.21) 

(B.22) 
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and 9 is the angle between _ and the z-axds. We shall work in the center of momentum 
frame of the two photons, and use the fact that for real photons, 


*1 • fcl = «2 * *2 = 0 . 


(B.23) 


After some algebraic manipulation, we obtain 


f{6) = /^sin 4 & 


in which 


+ 2 — 2/? 2 sin 2 0 


0 = P - 


[(1 — 0 COS 9) 2 (1 — 0 COS 0)(1 + 0 C06 8) (1 + /3 COS 

i i 


[1 — 0 cos 9 1 + 0 cos 8 


P-o 


(B.24) 


(B.20) cam be simplified by carrying out the integration over 8, so that 


L 


f{8) sin 0 = 4(2 - 0 2 ) + 2(0 2 + 1)(/? 2 - 1)^ In 

0 0 


1+0 


1-0 


(B.25) 


and hence 


a = 


e 4 a 2 n 4 pl 0 

° 0 x 2 tt 


(2tt) 10 • 2c 2 


1 ,11+01 


4(2 - 0*) + 2{0 l + 1)(0* - 1)- ln|^- 

Using the definition of a in (B.22), we then obtain 

\8n-4 1 - 2 


a = e 4 (2ir)®Q" 


16p 4 (2tt) 10 x 2c 2 


+ 2(^ 2 + 1)(0 2 — 1) ln|j-±' 


P—O n 
x — - x 2 tt x 


4(2 - 0 2 )0 


(B.26) 


(B.27a) 


wi&rz * 2w * [mi -^+2(^+ 


(2 -, 

1+0 


x In 


1-0 


(B.27b) 
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In the 4t natural units” in which one sets h = c = 1, this result can be written in the form 


' ■ = x 2 " + ^ + 1)(/?2 - *> * 5 Htt||] • (R28) 

In terms of the variables y = W = p~ 0 + P+o — 2p~o, "we can write 

<t = a 2 x ^ x 2 tt (1 + y)\/l _ V - (2 - y)y In - 1 (B.29) 

in which a = is the fine structure constant. This result is the same as the one obtained by 
Papageorgiu (ref. 2, eq. (15), p. 159). 77 — ► / + /~. For this case the interaction lagrangian 
can be written in the form 

Ant = — e rp(x) 4(i)^»(x), (B.30) 

in which i/>(x) denotes the lepton field operator. JL{x) = A^x) 7 ^ and 7 ^, p = 0, 1, 2, 3, are 
the Dirac 7 -matrices. (See ref. 9, Appendix 2 , p. 335-361.) tp(x) = yfi{x)i 0 , where ^(x) 
is the hermitian conjugate of xp(x'). The second order team in the 5 - matrix is defined by 
(B.15), with AntM defined by (B.30). The initial and final states can be denoted as 

It > = |fci»ci; fc 2 .C 2 >» (B.31a) 


/ > = p-,5_;p+,5+ >, 


(B.31b) 


in which we have already defined kj,tj,j = 1,2 as the photon momenta and polarization. 
p + ,p + are the momenta and spins of and l~ respectively. Following the notations of 
reference 9, appendix 2, we can write the 5-matrix element < /|5^|t > as 


< /|5< 2 )|t > = — e 2 a u(p_,5_) ^ v (P+> 5 +) 

+ u(p_,5_) fo-r?jr- — -5-7-T- 4 \ U (P+»S+) ^(P- + P+ ~ k l “ *= 2 )* 
k rz - m* + te J 


where 


(B.32) 
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k = p- - ki = ^2 - P+, k’ = p- - k2 = k\ - p+ , 


a = (2 tt) 4 (2 p-o • 2p +0 ■ 2fri 0 • 2 k 2 o)~ 1/2 n 2 , 

u(p_,5-) and V(p+, S+) are the spinor wave-functions associated with l~ and / + . This 
5-matrix element can also be represented diagrammatically by Figures. (B.la) and (B.lb). 
The total cross section is given by a formula similar to (B.19): 

S-f+ 

Performing the sum over the photon and lepton spins, we can write 

Y I </|s< J >|i>| 2 =£eW(p-+p + - + 

*1,£2 

S-,S+ 

x Tr[fi- + m) i{# + m) m) fo (-i)(/f + m) /ij 

+ (Jb 2 - m 2 + ie)~ 1 (k' 2 - m 2 + xe ) -1 Tr [(jf_ + m) i\ (k + m) / 2 (j4 - m) i\ (ft' + m) / 2 ] 
+ (fc 2 - m 2 4 - ic) - ^^ - m 2 + ie ) _1 Tr [(jf_ + m) / 2 (/f' + m ) /l(/4 “ m ) M# + m ) *l] 

+ (fc /2 - m 2 -I- ie )“ 2 Tr [(jf_ + m) / 2 (fc' + m ) i \ (?+ ~ m ) /l(^ + m ) * 2 ] }» (B.34) 

in which Tr denotes the trace operator. We use m instead of m/ to denote the lepton mass. 
From (B.34), it can be seen that the sum in (B.34) can be naturally divided into four terms: 

T\ = Tr [(j- + m) j\ {ft + m) / 2 (jf+ - m) + ™) /l] » (B.35a) 

T 2 = Tr [(] 6- + m) j\ ()i + m) /i(j4 - m) fo{jt r + m) / 2 ] , (B.35b) 

T 3 = Tr [(jrf_ + m) / 2 (^ + m) - rn) + m ) /l] * (B.35c) 

T 4 = Tr [(jrf_ + m) / 2 (^' +m) <i(j4 - m) /1 (#' + m) / 2 ] • (B.35c) 

Hence we can write 

Y I < /|s (2) |i > f = + P+ - *1 - *2) jUji 

‘l-<2- ' ' 

S-J+ 



2 d?p-d 3 p+ n2 n 1 

(2tt)6 2cT 0 . (B.33) 
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(B.36) 


x ^{(k 2 - m 2 + ie)~ 2 T\ + ( k 2 - m 2 + t€) _1 (fc' 2 - m 2 + ie) _1 (T 2 + T 3 ) 

£ 1, £ 2 

+ ( k 12 -m 2 + ie)~ 2 n}. 


After some straight forward though tedious mathematics, one arrives at the following: 

-Y r 1 =8fcf 0 (l-/Jcos9)(l + ^+ ^8 cos 9 + (Pooft) -8m 4 , (B.37) 

4 ri \ *io k 10 / 


£ i, £ 2 


- E r 2 = - E r 3 = 8 *V 2 (1 - «* 2 ») [l - ^ 2 (1 - ccs 2 0)] , (B.38) 


4 4 

£ 1, £ 2 £ 1, £ 2 


/ 2 2 \ 

- V T 4 = 8fcf 0 (l + /? cos 0) ( 1 + ^ 0 cos 0 - p 3 cos 3 0) -8m 4 , (B.39) 

4 ^ V *10 *10 ' 


£ 1> £ 2 


IT J 


where (5 = y—. Prom (B.33) and (B.36) we obtain 


- , -4n-4 _L 


" = ^ ^ 

x fl - Y, l (* 2 - rn 2 + «)- 2 T! + (* 2 - m 2 + is )- 1 (k n - m 2 + ie )' 1 (T 2 + T 3 ) 

+(fc ' 2 - m 2 + ie) - 2 T 4 ] jt^p-dV-i-^CP- + P+ - fcl - *2) 

+(Jfc 2 - m 2 + te) _1 (fc' 2 -m 2 + ie) _1 (T2 + T3) + (fc' 2 - m 2 + te) -2 T 4 ] j 

(B.40) 


x - sin 6 dd x 27r. 
2 


Now we use 


k 2 -m 2 = — 2fc 2 0 (l - 0 cos 9), 


(B.41) 


k 12 -m 2 = -2fc 2 0 (l + /3cos 0), 


(B.42) 
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together with (B.37), (B.38), and (B.39) to arrive at 

ri E (*- m ’ + *r* sin^^= 

•*0 ^ •'O € , .co 


C1,C2 


= -i/J 2 -4-8£.+2(2 + 3^ 

° *10 


m 2 1 


«1,«2 
1 + 0 


*10)^ 


In 


1-0 


4r?v 




(B.43) 


f - (fc 2 — m 2 + ie) 1 (fc /2 — m 2 + ie) 1 T'2 sin 6 dB 
«1>«2 

= f - (fc 2 — m 2 + tc)~ 1 (fc /2 — m 2 4- tc)~ 1 r3sing d.0 

Jo 4 £ l> e 2 

= 8(1 - 0 2 ) + i 0 1 - |(2 - /J 2 )(l - 0 2 )ln|{-^ 1 


(B.44) 


So from (B.43) and (B.44) we have 

f i ^2 [(* 2 ~ m2 + it)~ 2 T\ + (k 2 - rn 2 + tc) -1 (fc' 2 - m 2 + te) _1 (T2 + T3) 
■'O 


«1>*2 


+ (fc 72 - m 2 + te) -2 ^] sin 0 dB 
= -8(2 - 0 2 ) + 1(2 + 3(1 - /J 2 ) - (2 - /) 2 ) - (1 - 0 2 )]ln |i± 


1+0 
P 


— 8d + ») + ^(l + »-V)>“|^ + VP* 

where y = = 1 — /3 2 . Putting this into (B.40), we have 

-( 1 + y ) + |( 1 + y - y ) 


(B.45) 


1 1 


a = 


4-5: Wo x 8 


x In 


(2 tt ) 2 16fcj 0 2c 

Tv + JV~ 


2ir 


= x [- (1 + y)v/ ^ + 2 (‘ + " _ ^) to I ^ + Vr^ 

= ^ x 4 .[_ (1 + y)v ^ + 2( 1 + v-|) I n|-L + ^/p 


(B.46) 


in “natural” units. This result is the same as the one obtained by Papageorgiu (ref. 2, 
eq. (14), p. 159). 
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We now consider the process 77 — ► S° in which S° is a neutral scalar. Using pa to denote 
the momentum of S°, the cross section for this process can be written as 

12 n 1 


=/; = £ 

«1.«2 


<Pa 




-y ft X — X — , 

(2ir)^ 2c T 0 


(B.47) 


in which S denotes the 5-matrix. For the reverse decay process S° — + 77, the width T can 
be written in the form 


4/Eh *i» c i;*2,e2 

0 ^ £ 1, £ 2 ' 


Pa 


’<£*1 <^*2 „2 
(2-r) 6 


(B.48) 


From conservation of momentum, we can write 


<Pa 


ki,e;k 2 ,C2 > = <Ps 


ki,ei]k 2 ,€2 > 6*{p s - ki - k 2 ) (B.49) 


From (B.47) and (B.49) we now have 


= iE 


<Ps 

I 

e l. £ 2 

1 Q 2 fiT 0 
< (2tt) 3 2cT 0 (2tt) 4 

1 ft 3 




£(Pao ~ *10 - k 2 o) 


8c (27r) 


tE 


£ 1> £ 2 


<Pa 




fi(Pso ~ 2 A»i 


(B.50a) 


(B.50b) 


In (B.50), we ass ume that we are working in the next frame of S°. Likewise (B.48) can also 
be rewritten in the form 


r =^E 


£ l> £ 2 

ft 3 


< *i»€i;*2>«2 


Pa > 


'kl 


10 


x 4ir x 


fi 2 QT 0 


(2tt)9 


fc 10 51 

£ 1> £ 2 


< ki,ey,k 2 ,e 2 


2 (2tt) 6 (2tt) 4 

2 

1 

Pa > 


(B.51) 


From time-reversal invariance, we know that 


E 

< ki,ei;k2,e 2 

T 

Pa > 

2 

-E 

< Pa 

T 

*i»ei;fc2.€2 > 

«1,C2 




*1, *2 





(B.52) 
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Therefore from (B.50) and (B.51) we now have 

ff — — (27t) F X — j 6{pso 2fcio)* 

In the next frame of S°, 

2 

a = (27r) 2 r r” - 6(m 2 - 4fc 2 0 ) = ^-T6(m 2 - 5), 

«io m 3 

the square of the total momentum (k\ + J^) 2 . For 77 — + H° i 
particle, T can be written in the form given by (B.4) ^ (B.7). 


in which 5 is 
neutral Higgs 


(B.53) 


(B.54) 
which H° is a 
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Appendix C. Fermion Contribution to T(H° — ► 77). 

In this appendix, we derive the Fermion contribution to the decay width of the decay of 
a Higgs particle H° — ♦ 77 . For this case the interaction lagrangian can be written as (ref. 
13, eqs. (22.58), (22.78), pp. 676, 682) 


Ant( x ) ~ 4t(*> 4t(*)» 


(C.la) 


where 

4t( x ) = 4(*)ty(*).4t( x ) = h f i> f (x)7i(x)rf> f (x). (C.lb) 

■4> fix) is the fermion field operator, A^x) the photon field operator, and r?(x) the scalar 
Higgs field operator, qf denotes the charge of the fermion, and hj the coupling between the 
Higgs scalar and the fermion. The process H° — » 77 is third order in the interaction, so that 
the relevant term in the 5 -matrix is 


T^£ in t(x,)£ ill t(x2)Ant( x 3) d 4 *. d 4 X2 d^x 3 
The initial and final states can be denoted as 


(C. 2 ) 


t > = \ph >. and |/ >= |fci,ci;fc2,£2 >» 


(C.3) 


in which p/, denotes the momentum of the Higgs scalar, kj,ej,j =1,2, are the momenta and 
polarizations of the the photons. 

We use m j and to denote the masses of the fermion and Higgs scalar. The width for 
the process is given by 

|2 rf 3 *, d 3 k 2 

| < «l,Ci;K2,€2|*3'' , |P/i 


r - / E | < > f n 2 - (c.4) 
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The 5-matrix element can be represented by the diagrams: 


Employing standard techniques of field theory, we find 

< *l,€i;*2,€2| S (3) | Ph > = -q}hf(2pho • 2*10 ■ 2* 2 o)‘ 1/2 ^‘ 3/2 

x y^pijl^pi -*i) /2 -mj + t«] T r [(pi +mj) jti +mf) j 2 {t\ ~ #l~ #2 + ™f 

+ [(Pi -k 2 ) 2 -mj +ie] T r [(^i + m^) /2(^l~ fa + ™f) &+***/)] J 

x (pi - mj + ie) -1 [(pi - *1 - *2) 2 - mf + iej ^(fci + k 2 - Ph)- (C.5) 

We can separate the two terms on the right hand side of (C.5) and let = S^\ 

so that 

< *l,ei; *2, €2 1 s[ Z) \p h > = -q)hf{2p ho ■ 2* 10 • 2* 20 )~ 1/2 fr 3/ V(*l + k 2 -p h ) 

X y^PlfCPl-fcl ) 2 -^ 2 +*«] (Pi-mj + te ) _1 

x [(pi - *1 - *2) 2 ~ mj + tcj T{ 3) , (C.6a) 

where 

= T r + mf) /i(jrfi- ft + mf) ft- ft + my)j , (C.6b) 
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and 


< *i, ci; fe, e 2 | S®\-h > = -q}h f (2p ho • 2fc 10 • 2fc2o)‘ 1/2 + *2 ~ Pk) 

x J d^pi ^(pi - fci) 2 - ro 2 + * c ] [pi “ m / + * c ] 

x ^(pi - fcl - ^2) 2 - m / + ic ] ^2 

where 

T^ 3) = T r J(^i + m f ) fa+ m f) /l(^- ft" fa + m f)\ ( C - 8 ) 

The evaluation of the matrix elements (C.6) and (C.6) are quite similar. So we need only 
consider (C.6) in detail for illustration. By evaluating the trace in (C.6b), and using the fact 
that in the center of momentum frame of the two photons, 



ei • *1 = «1 • *2 = c 2 • *1 = c 2 • &2 = °> 

(C.9) 

and also 




?r 

►-**o 

II 

& 

n 

p 

(C.10) 

for real photons, we find 




t| 3) = 4m^(4pi • eipi t-t\- e 2 pj) + 8m/ei • *2Pl ‘ fc l 

+ 4rnf(-fci fc 2 + ™ 2 ) c l* € 2- ( C11 ) 

Now we use a standard technique of Feynman parameterization (ref. 14, Section 3.2, 
pp. 160-197). 

(p 2 - m 2 + it)* 1 [(pi - ti)“ - m* + «J [(PI - *:i - *2> 2 - ”> 2 + «] 

= 2 f\J^dy[(n-Q) 2 ^Y 3 ‘ C12a) 

where 

<?« = xk» + »(*? + *£), ( C12b ) 
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Pi = -Q z -m,f+ kfx + ( ki + k 2 ) z y + ie 
= -Q 2 - m 2 + 2ki ■ k$y + ie. 

Now (C.6) can be rewritten in the form 

< fci, «i; k2, C2 S 3 Ph > = —q 2 hj(2pf lo • 2k\Q * A:2o) _1 ^ 2 ^ _3 ^ 2 ^ 4 (^i + — Ph) 

x2 Jo dx Jo dy f d ~ Q) 2 + Pi] T ] 

Prom (C.ll) and (C.13), it is apparent that in order to evaluate (C.13) we need to compute 
the following integrals: 

I% u = J dx dy d*pi p^p\ [(pi - Q) 2 + p \ j , (C.14) 

dx dy d*pi pf [(pi - Q) 2 + p?] , (C.15) 

dx dy dtpi ^(pi - Q) 2 + p \ J (C.16) 

These integrals can be computed using the method of dimensional regularization in which 
one first computes the following integrals: 

I^{n) = J dx dy cPpi p^p\ [fo - Q) 2 + pf] , (C.17) 

7 i*( n ) = J dxdy (Ppi p%6 [(pi - Q) 2 + p \ J , (C.18) 

dx dy (PpiS^ipx - Q) 2 + pf] , (C.19) 

in which n is a real number, which in the final result are allowed to approach 4. Details of 
this process is given in Appendix D. 

From the results in Appendix D, we find 


(C.12c) 


i (3) . (C.13) 


1 [(pi-«) 2 +p?] 2p i l J 

pf 


/ 


d 4 


(Pi ~Q ) 2 + p\\ 


ii-eVw], 

(C.20) 

.IT 2 Q» 

2 p\ ’ 

(C. 21 ) 
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and 


I ** 1 


[<pi-<3) 2 +Pi] 


it 2 l 
3 = -Tp* 


(C-22) 


Define 


Jl = J d?Pi [<J>i-ki) 2 -m 2 f + ie] (p?-m}+ie) 
x [(pl ~ *1 “ *2) 2 ~ m / + tc] x t [ 3 \ 

By using (C-9)-(C-12), (C-20)-(C-23), we find 

Jl = —i7 r 2 J dx J dy -j j^my ^p 2 — Q 2 ) 

+ 8m jQ ■ k\ + 4my fci • &2 + m /)] C 1 “ c 2- 
Using the definitions of pi and <? in (C-ll) and (C-12), we can simplify (C-24) to 


(C-23) 


(C-24) 


Jl = —ITT 2 2m 


’ €l l *h *L 


1-1 4y(l - i - y) - 1 


dy 


y(l — x — y) — A + te 


(C-25) 


2 2 
m, mi 


771* Tr*# 

where A = sf- 

The integral in (C-25) can be done after some changes of integration variables and applying 
some techniques in complex analysis. The result is 

Jl = -i7r 2 2m f ei • £2 *A)> 

4 


(C-26) 


m'e s 

where for A = —*r > 7i 

mr “ 


J 0 = 2 - 2 (4A - 1) jarcsin » 


(C-27a) 


and for A < j, 


J 0 


= 2 + (4A-d{-y + 5 


In 


2 + 


-v/F 


— i7rln 


2 + 


2 ~ \[\ ~ A 


}■ 


(C-27b) 
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From (C-5), (C-6), (C-23) and (C-26) we have therefore 

< fci , Cl ; fc2» € 2|sj 3 ) | Ph > ~ ~9l ( 2 Pho ' 2^10 • 2 ^ 20 ) ^ ^ 3 / 2 Jl^ 4 (*l + *2 “ Pfc) 

= — g 2 hj (2ph a • 2fcio • 2*20 ) 1 ^ 2 0 (*i + *2 ~ Ph) 

x f-t7r 2 x 2m f x ci • C 2 J 0 ) • (0*28) 


It is straight forward to check that 

< k\, e\] k 2 ,e 2 ,\S 2 ^ \Ph > = < € li ^ 2 » € 2 |*^i ^|P/i- (C-29) 


Therefore using (C-3)-(C-8) and (C-28), 


r = — - f ^2 4 ?/ h ) ( 2p * 0 ‘ 2fcl0 ‘ 2fc2 °) 1 Q_3tf4 ( kl + k2 ~ Ph * 
T 0 4 J ei , £2 

A 9 , ,2,t.2 Wo d 3 fcld 3 fc 2 n 2 

x 7T 4 x 4m f x (C! • c 2 r I Jol x J2nf ' 


(C-30) 


Now 


Therefore 


( c l ’ C2 ^ 2 = 2 ’ 2fcl0 = 2 * 20 = Ph ° = m/l ’ 

ei,C2 


(C 31) 


r = igfe x 2 x — 3 ^ x 2 I Jol 2 J d?k\6 (*10 + *20 - Pho ) 


-4_,2 


= TeTW 1 Jo1 


(C-32) 


hj is related to the Fermi coupling constant Gp by h 2 — m z x 2\f2Gp (ref. 13, Section 
22.2, eqs. (22.58), (22.70), and (22.83), pp. 676, 679, 684) 

r = 16 xW ^ * 2 ' f2 ° F ^ ^ |J °' 2 

= ^ 4 S m * A2W2 ' <C33> 

We note that ^ is the charge of the fermion in units of the electron charge. Our result 
agrees with that in the literature. (See ref. 12, eqs. (10) and (11), p. 95.) We note that the 
sign of the imaginary part of J Q in (C-27b) is opposite to that of reference 12, equation (11). 
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However, since only |J 0 | 2 enters into quantities of physical interest, such as T and <r, therefore 
this difference in sign of the imaginary part of J Q is not significant. 
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Appendix D. Evaluation of certain integrals. 

In this appendix, we outline the procedures involved in evaluating the integrals in (C-14)- 
(C19). 

The following integrals can be evaluated by standard methods of calculus. 


“ “ mdu 


r- 

Jo (« 


(Dl) 


x 2 + M 2 ) 1 

provided m is even, m > 0, l is an integer > 0, t > y + 1, and the coefficients C m ^ are 
defined by 

m 


C m,l ~ y 


(£ — y ~ 1)! /m — 1 m — 3 3 I 


T ZH ( 

-1)! i, 


2 2 ) 


t—m—2 


X /^(cosfl) 2 *" 

Jo 

If t is a half-integer, (D l) stiD applies with 

x / ? (cos fl) 2 *" 

^0 


(D-2) 


£— m— 2 


(D-3) 


Using methods of complex analysis, we can show that the same formula (D-l) applies if 
M is replaced by iM in (D-l). 

f <Tp(p 2 + m 2 )-“ = r° ip „ j°° lp|"- 2 <i|p| (p 2 - |p| 2 + M 2 )-«n™ 2 , (D-4) 

where denotes the surface area of the n-dimensional unit sphere. Now using (D-l), we 
find 


J <Pp{j> 2 + M 2 )~ a = dp 0 (-l) a C n . 2 ,a(-Po - M 2 )- a+a T 


n-1 


= 2(-l)^C„_ 2 , a C 0ia _ 2 - 1 f!< 0 ! 2 M- 2a+ ". 


(D-5) 
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From (D-5), it is straight forward to compute 


J <Pp(p 2 + 2pQ + M 2 ') “ = j d n p[(j> + Q) 2 + M 2 - C? 2 ] 

= 2(-l)T 1 C n -2 A C o,a-^ n i-2 (" 2 - &) 


— a-H 


(D-6) 


Now we can evaluate 


/^(p^p-Q + M 2 )'^^ ^/*"p(p 2 + 2p.<? + M 2 ) 


-a+1 


= (-i) 2 *" 


5+1 2 (o - 1 - 5) 


a — 1 


<^ 21 - 2 , 0-1 


X C^aj. X n™ 2 (a / 2 - Q 2 ) '“ +5 <?„. (D-7) 


Using similar techniques we can evaluate 


J cTp PtiPv (p 2 + 2p • <3 + M 2 ) a = (-l) S ^(a-2-^)(a-l)- 1 (a-2) 1 

x Cn— 2,a— 2 ^ a _^ (m 2 - g 2 )"° +1+5 

X 3 M „+(M 2 + Q 2 ) -1 (-a + l + £) 

x (-2g M Q„)] , (D-8) 


and 


/ <fp p 2 (p 2 + 2p • Q + M 2 ) “ = (-1) 3 ? (a - 2 - (a - l) _I (a - 2)- 1 

X (m 2 - < J 2 )'“ +1+3 


n 


= 2 (“- 1 -f) K-«T‘« 


, (D-9) 


in which q^, is the metric tensor 


9oo = “911 = -922 = -933 = 1 - 


(D-10) 


and all g^u with p ^ v are 0. 
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We also note the following: 


f°° u m -i a r°° u m 

Jo (u 2 + M 2 )“ U_ o-l 0 M 2 /o (uS + A^JO - 1 

Therefore, by using (Dl), we find 


du 


Cm, a — — 1 2 — ^ ( a — 1) 


We can now use these results to evaluate 


I <r Pl (4 pfrf - P ?<r) [(pi - Q) 2 +p?] 0 

= (-l) 2 ? 1 (a-2-2) (a- l) -1 (a - 2)- 1 C n _ 2 , a -2C 0 a _n+3^°2 2 


<T (4 - n) (p? + Q 2 ) (p?) ° 2 - 4(4 - n)Q ,i Q v (pj) 
Using (D-12) in (D-13), we find 

J <Ppi ( 4 PiPi - PiP^) [(pi ~ Q) 2 +pi] 

= (-1) 2 * 2 (a - 1 - 5 ) _1 c n _ 2 , a c 0ja _^ni 0 2 2 


-a+3 


-a+5 


<7^(4 - n) (p? + Q 2 ) (p?) “ 2 - 4(4 - n)<W (p?) 

Setting a = 3, and taking the limit as n — ♦ 4, we find 

| <f w (4pf p f - P ?<r) [(p, - q) 2 + P f] 

= (-1)5 x 2 nfi [<r (p? - Q 2 ) + w-] c 2 , 3 c 0 j 

= -iyi[r ( P f - q 2 ) + . 

In similar fashion we find 

J [(PI - Qf + Pi] ~ d = (-l) 2 ^ 2 (a - 1 - ^) (o - 

= (-i ) 2 ^ 2 2 c„_ 2 ,„ c n™, (pf)"“ +? (-<?*) 


(Dll) 


(D-12) 


(D-13) 


(DM) 


(D 15) 


(D-16) 
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Therefore, 


J d*Pl Pi [(Pi “ Q ) 2 + Pl] — - * 2< ^2,3 Cq,| ^2 ^ p 2 Q * 1 




, *2 QP 
2 pf 


Finally, 


J <^PlPi [(Pi - Q) 2 +Pl] — 2 ( _1 ) 3 / 2 C2,3 C 0 3 


. tt 2 1 

= -t “o" “5- 
2 Pi 


(D-17) 


(D18) 
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Appendix E. Stopping Power 

Consider the reaction 

Z\Z 2 — » Z1Z2X, (E-l) 

in which X represents one or more particles produced in the process. Let Z\ be an incident 
particle, and Z2 represent a fixed target, whose density is p (number of nuclei per unit 
volume). Let a denote the cross section for the process (E-l), and E x the energy of the 
system X. If we disregard the effect due to recoil of Z2, then by the conservation of energy, 
the energy loss of Z\ is equal to E x . Consider a slab of the target Z z of cross-sectional area 
A and thickness Ax. 


Figure E.l 


The number of Z2 nuclei in this slab is pAAx. The cross section for an incident particle 
Z\ to collide with a Z2, producing X is given by 

A a = pAAx dE x , (E- 2 ) 

ihfx 

where we assume the energy of the produced system X to be between E x and E x + dE x . 
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Therefore the probability for this process is 


P(E x )dE x = — = pAx dE x , (E-3) 

in which P(E X ) represents the probability density for the process. Therefore the total energy 
loss by the incident particle Z\ per unit length is given by 



(E-4a) 

(E-4b) 


The — sign in (E-4) signifies the fact that energy is lost by Z\ in the process, so that the 
change in its energy dE is negative. (See ref. 15, eq (6-4), page 741. 

For two-photon processes of this kind that we have considered 


E x = w i + w 2 , (E-5) 

in which again we use the “natural units” for which h = 1. The cross section is given by 
(2-2a). By switching the variables of integration from w\, W 2 to w\, E x , and using the fact 

dw\ dw 2 — dw\ dE x , (E-6) 

which can be obtained from (E-5), (E-4b) can be written in the form 
dE f dw i f dE x 

—fa = P J — J - E ' x - Wl Ex F ( Wl ' Ex ~ Wl ^ x <T 'rr( w i’ Ex ~ ^i) ( E ' 7a ) 

= p /— (tiJi +iu 2 ) F{wi, w 2 )(7 r y(wi, w 2 ) (E-7b) 

J Wi J W2 

(E-7b) is the same as (2-7). 
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Figure Captions 


Figure 1. 

Figure 2. 
Figure 3. 
Figure 4. 

Figure 5. 

Figure 6. 
Figure A.l: 
Figure A. 2: 

Figure A.3: 
Figure B.la-b: 
Figure B.2: 
Figure C.la-b: 
Figure E.l: 


A Feynman diagram for the processes Z}Z 2 — ► ZiZ 2 l + J , ZjZ 2 — ► ZiZ 2 s + s 
and Z X Z 2 — Z 1 Z 2 V + V r ". 

A Feynman diagram for the process Z\Z2 — ► ZiZ 2 H*\ 

Plots of W 2 -jfa. 

a-d: 0 'yt(W 2 ) for the reactions 77 — * l + l~, 77 — * s + s“, 77 — ♦ V + V~ anc 
77 —* H® . 

Plots of P (b) for the reaction 208 P6 208 P6 — * 208 P6 208 P6e + e - at different 
energies. 

Plots of the total cross section for the process Pb 2 ®* Pb — > 208 P6 208 P6P°. 

Electromagnetic fields generated by a charge q moving along the x-axis. 

Protons emitted by two colliding nuclei, viewed along direction of motion of 
the nuclei in their center of momentum frame. 

Cro6S-sectional view of the collision of two nuclei. 

Second order Feynman diagrams for the process 77 — * s + s“. 

First order Feynman diagram for the process 77 — ♦ s + s“. 

Feynman diagrams representing fermion contribution to the process — ► 77. 

A beam of particles Z\ incident on a fixed target Z 2 . 
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